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Abstract 

QJ , In the unconstrained MSSM, we reanalyze the constraints on the phases of supersymmetric flavour 
conserving couplings that follow from the electron and neutron electric dipole moments (EDM). We 
find that the constraints become weak if at least one exchanged superpartner mass is > 0{1 TeV) 
^ ■ or if we accept large cancellations among different contributions. However, such cancellations have 

» . no evident underlying symmetry principle. For light super partners, models with small phases look 
like the easiest solution to the experimental EDM constraints. This conclusion becomes stronger 
the larger is the value of tan/3. We discuss also the dependence of ex, ^niB and b ^ sj decay on 
those phases. 
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1 Introduction 

In the Minimal Supersymmetric Standard Model (MSSM) there are new potential sources of 
the CP non-conservation effects. One can distinguish two categories of such sources. One is 
independent of the physics of flavour non-conservation in the neutral current sector and the 
other is closely related to it. To the first category belong the phases of the parameters fi, 
gaugino masses Mj, trilinear scalar couplings Ai and m^25 which can in principle be arbitrary. 
They can be present even if the sfermion sector is flavour conserving. Not all of them are 
physically independent. 

The other potential phases may appear in flavour off-diagonal sfermion mass matrix 
elements Am?- and in flavour off-diagonal LR mixing parameters Aij. These potential new 
sources of CP violation are, therefore, closely linked to the physics of flavour and, for instance, 
vanish in the limit of flavour diagonal (in the basis where quarks are diagonal) sfermion 
mass matrices. It is, therefore, quite likely that the two categories of the potential CP 
violation in the MSSM are controlled by different physical mechanisms. They should be 
clearly distinguished and discussed independently. 

Experimental constraints on the "flavour-conserving" phases come mainly from the elec- 
tric dipole moments of electron [0 and neutron 0Q: 

E'fP < 4.3 • lO^^^e ■ cm 
E'fP < 6.3 • 10-26e . cm 



n 



The common belief was that the constraints from the electron and neutron electric dipole 
moments are strong [^ ^ and the new phases must be very small. More recent calcula- 
tions performed in the framework of the minimal supergravity model |^, |^, |^ indicated the 
possibility of cancellations between contributions proportional to the phase of fi and those 
proportional to the phase of A and, therefore, of weaker limits on the phases in some non- 
negligible range of parameter space. The possibility of even more cancellations have been 
reported in ref. [§ in non-minimal models. For instance, for the electron dipole moment, the 
coefficient of the /i phase has been found to vanish for some values of parameters. Since the 
constraints on the supersymmetric sources of CP violating phases are of considerable theo- 
retical and phenomenological interest, in this paper we reanalyze the electric dipole moments 
with the emphasis on complete understanding of the mechanism of the cancellations. 

The new flavour-conserving phases in the MSSM, beyond the Oqcd present in the SM, 
may appear in the bilinear term in the superpotential and in the soft breaking terms: gaugino 



^Status of the new experimental limit on the neutron EDM is still under discussion 



masses and bi- and trilinear scalar couplings - see eqs. (X]3,^]|). We define them as: 



/^ „i<*. _ Mi ^,. _ _ Ai .. _ m 



12 



|/i| |M(| |A/| \m 
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Not all of those phases are physical. In the absence of terms ([A.3| , |A.6|) the MSSM Lagrangian 



has two global f/(l) symmetries, an R symmetry and the Peccei-Quinn symmetry p!0[ 



Terms ( |A.3| , [A.6[ ) may be treated as spurions breaking those symmetries, with appropriate 



charge assignments. Physics observables depend only on the phases of parameter combina- 
tions neutral under both f/(l)'s transformation. Such combinations are: 

Mi/i(m?2)* Ani{m\^Y A\Mi (2) 

Not all of them are independent. The two f/(l) symmetries may be used to get rid of two 
phases. We follow the common choice and keep m^2 ^^^1 i^ order to have real tree level Higgs 
field VEV's and tan/30. The second re-phasing may be used e.g. to make one of the gaugino 
mass terms real - we choose it to be the gluino mass term. 

A particularly simple picture is obtained assuming universal gaugino masses and universal 
trilinear couplings A/ = A at the GUT scale. In this case ?7(l)'s invariant parameter 
combinations (0) contain at that scale only two independent phases. Defining 0^ = 01 = 
02 = 03, 0A = 0Ai we can write them down as 

h = <\>A- 03, (^2 = 03 + 0/. - arg("^?2) (3) 

The re-phasing freedom may be used in this case to make all Mj simultaneously real. The 
RGE for Mi at one loop does not introduce phases once they are set to zero at GUT scale. 
The second f/(l) rotation can be used again to remove phase from m\2 already at Mz scale. 
Then only yU and Aj parameters remain complex at electroweak scale. Phases of various Ai 
parameters are not independent and can be calculated from the RGE equations. 

In most of the calculations in the next sections we keep in general /i, Mi,M2 and Ai 
complex. As one can see, all the physical results depend explicitly only on the phases of 
parameter combinations (§), as follows from the general considerations above. 

In Section 2 we discuss in detail the electron electric dipole moment. First, we present 
the results of an exact calculation, which is convenient for numerical codes. For a better 
qualitative understanding, we also perform the calculation in the mass insertion approxi- 
mation. The results of the two methods can be compared by appropriately expanding the 
exact results for some special configurations of the selectron and gaugino masses. After those 



^Loop corrections to the effective potential induce phases in VEV's even if they were absent at the tree 
level. Rotating them away reintroduces a phase into the m\2 parameter. 



technical preliminaries we discuss in Section 2 the magnitude of various contributions to the 
electric dipole moment and investigate the pattern of possible cancellations. The first obser- 
vation we want to emphasize is that, even without any cancellations, there are interesting 
regions in the parameter space where the phases are weakly constrained. Secondly, we do 
not find any symmetry principle that would guarantee cancellations in the regions where the 
phases are constrained. Such cancellations are, nevertheless, possible by proper tuning of 
the /i and A phases to the values of the soft masses. 

In Section 3 we analyze the neutron electric dipole moment with similar conclusions. In 
Section 4 we discuss the role of the n phase in the ex measurement and in the b —>■ sj decay. 
The necessary conventions, Feynman rules and integrals are collected in the Appendices. 

2 Electric dipole moment of the electron 

2.1 Mass eigenstate vs. mass insertion calculation 

The electric dipole moments (EDM) of leptons and quarks, defined as the coefficient E of 
the operator 

Ce = -^E7/;c^^,75V^F^^ (4) 

can be generated in the MSSM already at 1-loop level, assuming that supersymmetric pa- 
rameters are complex. 

In the mass eigenstate basis for all particles, two diagrams contribute to the electron 
electric dipole moment. They are shown in Fig. |I] (summation over all charginos, neutralinos, 
sleptons and sneutrinos in the loops is understood). The result for the lepton electric dipole 





Figure 1: Diagrams contributing to lepton EDM. 



moment reads: 



em^ ^ ^ 



o^ j=l K=l 



167r2 



^ ^ m^ Jm ((\/;^ 



' ILNfL^yiLNyR'') Cl2{ml m\) 



(5) 

where (V^i>c)L, (VJi^c)/?, (VJ^Af)i' {^iln)r are, respectively, the left- and right- electron- 
sneutrino-chargino and electron-selectron-neutralino vertices and Cii,Ci2 are the loop in- 
tegrals. Explicit form of the vertices and integrals can be found in [Appendix A . 

The eq. (^ is completely general, but as we discussed already in the Introduction, in 
the rest of this paper we assume no flavour mixing in the slepton sector. Therefore, in the 
formulae below we skip the slepton flavour indices. 

We present now the calculation of the electron EDM in the mass insertion approximation, 
for easier understanding of cancellations of various contributions, and then compare the 
two results. We use the "generalized mass insertion approximation", i.e. we treat as mass 
insertions both the L-R mixing terms in the squark mass mixing matrices and the off-diagonal 
terms in the chargino and neutralino mass matrices. Therefore we assume that the diagonal 
entries in the latter: |/i|, |Mi|, IM2I are sufficiently larger than the off-diagonal entries, which 
are of the order of Mz- 

There are four diagrams with wino and charged Higgsino exchange, shown in Fig. 0. Their 






Figure 2: Chargino contribution to lepton EDM in mass insertion expansion. 



contribution to the electron EDM is {Ef. = E} for electron): 

{Ee)c - 



2^^^"^^Im(M.,)tan/3M^I^'-^)-^-(l^^l^'-^) 



(47r) 



/ip - IM2 



(6) 



Neutral wino, bino and neutral Higgsino contributions can be split into two classes: with 
mass insertion on the fermion or on the sfermion line. Diagrams belonging to the first class 
are shown in Fig. ^. Their contribution has a structure very similar to that given by eq. (^: 
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Figure 3: Neutralino contribution to lepton EDM in mass insertion expansion: mass insertion 
on the fermion line. 



{Ee)Nf = - Im(M2/i)tan/3- 



2(47r)2 ' "^' ^ |/i|2-|M2|2 

+ ^ Im(Mi/x)tan/?- '^^ £^' i'"! ^ i^v eh im 
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2(47r)2 ' '^^ " |/Up-|Mi 
~ (A \2 Im Mi/i tan/? ^ — — 7 

where we denote by rriE, mE'^ and mj^ the masses of left- and right- selectron and electron 
sneutrino, respectively. Finally, diagrams with mass insertions on the selectron line are shown 
in Fig. ^. Only the first two with bino line in the loop give sizeable contributions. The other 
two with neutral Higgsino exchange are suppressed by the additional factor O (ml/M^) and 
thus completely negligible. The result is: 

{Ee)Ns = — — -^Im[Mi(/itan/3-F A^^' 



(47r)2 ^ m^ — m 



E'^ 
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+ terms suppressed by C , 2 



M^^ 
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Figure 4: Neutralino contribution to lepton EDM in mass insertion expansion: mass insertion 
on the scalar line 



Eqs. (^§) have a simple structure: they are linear in the CP-invariants (Q), with coeffi- 
cients that are functions of the real mass parameters. Thus, the possibility of cancellations 
depends primarily on the relative amplitudes and signs of those coefficients. An immediate 
conclusion following from (|^-§) is that limits on the Mj/x phases are inversely proportional 
to tan/?. Therefore, in the next Section, we discuss limits on sin 0^ tan /5 rather than on the 
/i phase itself. 

The approximate formulae (|^-§|) work very well already for relatively small |/i|, |Mi| and 
IM2I values, not much above the Mz scale. Figure ^ shows the ratio of the electron EDM 
calculated in the mass insertion approximation to the exact 1-loop result given by eq. (^. 
The accuracy of the mass insertion expansion may become reasonable already for |yu| > 150 
GeV (depending on \M2/fi\ ratio) and becomes very good for \fi\ > 200 — 250 GeV. 

Formulae (P-||) can also be obtained directly from the exact expression (|^) using the 
expansion of sfermion and supersymmetric fermion mass matrices described in [Appendix B| 
(this gives a very useful cross-check of the correctness of the calculations). One should 
note that even though the expansion ( p.27|) does not work for degenerate IM2I and |/i|, 
the expression (P) has already a well defined limit for IM2I = |/x|. The same holds for the 
I Ml I = |ju| and the degenerate sfermion masses in eqs. (|^,|). 
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Figure 5: Ratio of the electron EDM calculated in the mass insertion approximation to the 



exact 1-loop result. Thinner lines: 



0, thicker lines: 0a = 0. Degenerate left and right 



slepton masses and Mi/a\ = M2/a2 assumed. 



2.2 Limits on /i and A^ phases 

It is useful to consider two classes of models: one with Mi/ai = M2/a2 for gaugino masse^, 

that is universal gaugino masses at the GUT scale (the universal phase can be set to zero by 

convention), and the other with non-universal gaugino masses and arbitrary relative phase 

between Mi and M2. In the universal case we choose /i and A^. phases as the independent 

ones, in the second case the Mi,M2 phases are the additional free parameters. Constraints 

from the RGE running and proper electroweak breaking are insignificant at this point, as we 

have enough additional parameters to satisfy them for any chosen set of low energy values 

for /i. Ml, M2, rriE and Ag 

We shall begin our discussion by presenting the magnitude of each contribution (^, (|^) 

and of the /i and Ae terms in eq. (H), separately. For the Mi/ai = M2/a2 case a sample 

of results is shown in Fig. ^. We identify there the parameter region where at least one of 

the terms is such that for sin c/)^ tan /3 fixed at some assumed value, its contribution to E^. is 

larger than £^|^p. Barring potential cancellations, the fixed value of sin 0^ tan/? is then the 

limit on this phase in the identified parameter region. In the left (right) plot of Fig. ^ we 

show the regions of masses (below the plotted surface) where the limits on | sin0^| tan/5 are 

stronger then 0.2 (0.05), respectively. We see that this region extends up to 900 (400) GeV 

■^Here, by ai, a2 we understand the gauge couplings in GUT normalization, i.e. without a 5/3 factor for 
the C/(l) couphng. 



in tue (physical left selectron mass) for small values of |/i| and |Mi|, IM2I (up to 200-300 
GeV, say) and for larger values of |yu| and/or |Mi|, IM2I it gradually shrinks to 100 (0) GeV 
in mE for |/i| ~ IM2I ~ 1 TeV. We assume left and right slepton mass parameters equal, 
Ml = Me, so that the physical masses of the left and right selectron differ by D-terms only. 
The regions below the plotted surfaces are the regions of interest for potential cancellations. 
We observe, however, that even without cancellations, there are interesting regions of small 
\fi\ and IM2I and rriE > 0{lTeV) or small rriE and |yu| ~ IM2I > (9(500 TeV) where the 
phase of /i is weakly constrained. One should also note that for very large |/i| and the other 
masses fixed the limits on the fi phase get stronger again. This is due to the term (^, which 
does not decouple for large |/i|. The limits on {\Ae sin (p j^ J /rriE are typically significantly 
weaker. 
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Figure 6: Regions (below the dark surface) for which generic limits on |sini 

stronger then, respecively, 0.2 (left plot) and 0.05 (right plot). Degenerate left and right 

slepton masses and Mi/ai = M2/a2 assumed. 



In Fig. we show again the limits on /i phase, this time as a two-dimensional plot in the 
[jnE-, 1^2!) plane, assuming Mi/ai = M^ja^ and ^a^ = 0. The limits plotted in Fig. |^ are 
given by the sum of all terms (|^-§|), not by the largest of them like in Fig. ^ 

In Fig. ^ we show similar limits on the Ae phase on {mE,\Mi\) plane, assuming vanishing 

/i phase. The limits on the A^. parameter phase are significantly weaker and decrease more 

quickly with increasing particle massesQ. They are almost independent of M2 and /i, as can 

be also seen immediately from eq. (H). 

'^One should note that the off diagonal entry in the s lepton mass matrices describing LR-mixing is propor- 
tional to Ae and to the electron mass me (see eq. ( A. 12 )). Therefore, even if imaginary part oi Ae parameter 



is weakly constrained, ImAe/mE ^ 1, the full LR-mixing in the selectron sector can have only very small 
imaginary part, of the order oimemE- For the electron, Gabbiani ct. al |ll[] give the limit IvaAe/mE < 10^^, 
but their definition of Ae contains me in it. After extracting it, their limit is similar to ours. 
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Figure 7: Limits on |sin0^|tan/3 given by the electron EDM measurements, sin^^e = 
and Mi/ai = M2/a2 assumed. 
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Figure 8: Limits on lAe/rriE sin (p^J given by the electron EDM measurements (|/i| 
GeV, sin0^ = and Mi/ai = M2/a2 assumed). 
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The magnitude and signs of individual contributions as a function of rriE are illustrated 
in Fig. ^ We plot there the coefficients of /i and A^ phases obtained from the exact 1-loop 
result and normalized by dividing them by the experimental limit on the electron EDM. 
Their shape depends mostly on the r?7,£;/|yu| ratio, much less on the \fi/M2\ ratio and scales 
like l/m"^. We see that either the chargino contribution to the term proportional to the 
II phase dominates (for small \fi\), or, if they become comparable (possible only for larger 
values of |/i| > 700 GeV), the chargino and the dominant neutralino contribution, given by 
eq. (§), to the /i phase coefficient are of the same sign. The neutralino contribution given by 
eq. (|^ has opposite sign than that of eq. (||), but their sum is positive. The only exception 
is the case of \fi\, \Mi\ < 100 GeV, where both neutralino contributions are much smaller 
than the chargino one. Thus, the full coefficient of the fi phase cannot vanish and the only 
possible cancellations are between the A^ and /i phases. Since the Ae phase coefficient is in 
the interesting region much smaller such cancellations always require large Ae in the selectron 



sector, Ae/rriE S> 1. This is shown in Fig. |T0], where we assume "maximal" CP violation 
0M = <^Ae = 7r/2. 

Better understanding of the /i-Ag cancellation can be achieved after some approximations. 
For light supersymmetric fermions, significantly lighter than sleptons, chargino exchanges 
dominate (Fig. 0), whereas in the opposite limit the biggest contribution is given by the 
diagrams with bino exchanges (Fig. |). Eqs. (H-|D can be greatly simplified in both cases, 
giving for degenerate slepton masses rriE ~ rriE'^ ~ "^jj: 

1) |Mi,2|, Ia^I ^ fTT-E- 

_ eg'^rrie Im(M2/i) tan/3 |/ip eg'^melm{MiA*) 

" "^ (47r)2 m|(|/i2| - |M2|2) °^ jMaP 2(47r)2 m% ^ ' 



|Mi,2|, l/il > mE. 
Ee 



egrrie Im(M2yU.) tan/5 eg' rrif. Im(Mi/i) tan/5 
4(47r)2 |/i|2|M2|2 4(47r)2 |^|2|Mi 



2 



/2 



efyf' me Im [Mi(/itan/5 + y4g)] / m 



5 + 21og— ^ (10) 



2(47r)2 |Mi|4 \ ° |Mi|2^ 

The behaviour of the lepton EDM is different in both limits. For heavy sleptons, |Mi^2|, |/^| "^^ 
iriE the coefficient of the /i phase decreases with the increasing slepton mass as l/m\. The 
coefficient of the A^, phase decreases faster, as l/rn^E- Therefore, in this limit the exact 
cancellation between A^, and /i phases requires large A^. value, growing with increasing niE- 
However, because all contributions simultaneously decrease with increasing m^;, partial can- 
cellation between /i and A^ phase is already sufficient to push the electron EDM below the 



experimental limits, what may be observed in Fig. |10| as a widening of the allowed regions 
for large uie- 
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Figure 9: Relative signs and amplitudes of various contributions to the electron EDM, 
normalized to (divided by) the experimental limit. Solid, dashed, dotted lines: coefficients 
of sin 0^ tan/? given by chargino (eq. (||)) and neutralino contributions (eqs. (^ and (|)) 
respectively. Dotted-dashed line: coefficient of |y4e sin0Ael/''^-E;- Degenerate left and right 
slepton masses and Mi/cti = M^ja^ assumed. 
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Figure 10: Regions of rriE — Mi plane allowed by the electron EDM measurement assuming 

"maximal" CP violation 0^ = 0^^ = 7r/2 and the corresponding values of Ae/rriE parameter 

(marked on the plots), necessary for the fj, — A cancellation. 



For sufficiently small slepton masses the full cancellation between fi and Ag terms occurs 
approximately for 



sin0^Jy4e| = sin0^|/i| tan/? 1 + 



iMil 



l/^r 5 + 21og 



m| 



:ii) 



where we assumed Mi/ai = M2/a2 and redefined the Mi phase to zero. Since this result 
is valid for |/i|, |Mi|, IM2I ^ rriE we see that for comparable 0^ and (pAe the cancellation 
is again possible only for large Ae/rriE ^ 1. For large |/i| ^ \Mi\, when one can neglect 
the second term in the parenthesis in eq. (0), the A^, giving maximal cancellation is almost 



independent of |Mi|, what can be observed in the right plot of Fig. |T0|. The allowed regions 
also widen with increasing |Mi|, but slower then for large niE because the /i and Ae phases 
are in this case suppressed by lower powers of |Mi|: l/|Mi| and l/|Mip respectively, instead 

of l/m'^E ^^'^ ^/"^E- 

It is worthwhile to note that in the most interesting region of light SUSY masses, where 
the limits on phases are strongest, the cancellation between (fixed) fi and A^ phases may 
occur only for very precisely correlated mass parameters, i.e. it requires strong fine tuning 
between |/i|, |Mi|, IM2I and \Ae\. 

In Fig. |Tl| we plot the allowed regions of the (pA^ — 0^ plane for chosen fixed values of 
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mass parameters. For light SUSY masses they are very narrow. This means that for fixed 
hght mass parameters one needs strong fine tuning between the phases in order to fulfill 
experimental limits. 
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Figure 11: Illustration of the fine tuning between the phases for fixed mass parameters 
(listed in the plots), tan/5 = 2 and two values of Ae/rriE = 5, 10 (marked near the respective 
allowed regions). Degenerate left and right slepton masses and Mi/ai = M2/a2 assumed. 



We shall discuss now the general case, with non-universal gaugino masses. The results 
for the magnitude of individual terms remain qualitatively similar. This is shown in Fig. |12 
for Ml = 100 GeV. The region of strong constraints on the /i phase shrinks in mE with 
increasing Mi. One can also see again some subtle effects like the expansion in tue of this 
region, for fixed Mi and M2 and increasing \fi\. This behaviour can be easily understood 
from the analytic results of eq. (H), where one can identify the term increasing with |/i|. 
The magnitude of individual contributions as a function of rriE has very similar behaviour 
as in the universal case - again, for small \fi\ chargino contribution dominates for all values 
of |Mi| and IM2I. The only possible cancellations for this |/i| range are between /z and Ae 
phases. For larger values of |/i| > 700 GeV the magnitude of individual terms may become 
comparable. For instance the /i phase coefficients in Ec and E^g terms (eqs. (||,||)) becomes 
comparable for \fi\/mE fixed by the ratio IM1/M2I. With arbitrary relative phase of Mi and 
M2 it is possible to cancel the terms proportional to the /i phase. To study this possibility it 
is more convenient to plot the contributions proportional to Im(yuMi) and Im(/iM2) (Fig. |T3|). 
They are comparable for r7i£;/|/i| ~ 1/5 — 1/3, depending on IM1/M2I ratio. It is clear that 
choosing 0i and 02 phases such that sin(0^ + ^2) and sin(0^ + 0i) have opposite signs, e.g. 
01 — 02 ~ TT, would give cancellation at these points. 
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10D 

Figure 12: Regions for which generic hmits on | sin0^| tan/3 are stronger then, respectively, 
0.2 (left plots) and 0.05 (right plots). The plots are done for Mi = 100 GeV, degenerate left 
and right slepton masses are assumed. 
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Figure 13: Relative signs and amplitudes of various contributions to the electron EDM, 
normalized to (divided by) the experimental limit. Solid, dashed and dotted lines: coefficients 

of sin(0^ + 02) tan/3, sin(0^+0i) tan/3 and |Ae|/m£;sin(0^^ — 0i) respectively. |Mi| = IM2I = 
100 GeV and degenerate left and right slepton masses assumed. 



14 



"M 



7r/2 (one can 



To give a more specific example, lets assume 02 
always achieve 

applies for |0i — ip2\ = iv^ii ^ /i/z,. rui y^iy ngm i-iwii cluu //tij, |(/^i 
order to suppress the term proportional to Im(/iMi). For heavier Mi and m.^; to avoid limits 
on 0^ one needs |0i| ~ vr - otherwise the Im(/iMi) term is suppressed too strongly. This 



and maximal 
2 = by field redefinition). In such a case, this new possibility of cancellation 

- 02| = 101 1 > 7r/2. For very light |Mi| and tjie, |0i| ~ vr/2 is required, in 



behaviour is illustrated in Fig. 14 
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Figure 14: Allowed regions for Mi phase as a function of IM1/M2I for some choices of mass 
parameters, tan/5 = 2 and 0^ = tt/2. Thicker lines: |/i| = 200 GeV, thinner lines: |/i| = 1000 
GeV. IM2I = 200 or 1000 GeV (marked on the plots), 02 = (f)Ae = and degenerate left and 
right slepton masses assumed. 



3 EDM of the neutron 



3.1 Formulae for the neutron EDM 

The structure of the neutron EDM is more complicated then in the electron case. It can be 
approximately calculated as the sum of the electric dipole moments of the constituent d and 
u quarks plus additional contributions coming from the chromoelectric dipole moments of 
quarks and gluons. The chromoelectric dipole moment Cg of a quark is defined as: 



J^c = --Cgqa^.^.T'^qG^' 



The gluonic dipole moment Cg is defined as: 



9 — ~'^^aJabc^ap^v ^\a^ 



fiuXa 



(12) 



(13) 
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As an example, in Fig. 15 we list the diagrams contributing to the (i-quark electric dipole 
moment. 
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d' d' N^ d' d' ~9 

Figure 15: Diagrams contributing to (i-quark EDM. 




Exact calculation of the neutron EDM requires the full knowledge of its wave function. We 



use the "naive" chiral quark model approximation |12], which gives the following expression: 



47r 



a) + ^^^a 



(14) 



3 ' ' 47r ' 'Air 

where rji and Ax are the QCD correction factors and chiral symmetry breaking scale, re- 
spectively, rje ^ 1.53, rjc ^ rig ^ 3.4 [|l3l. Ax = 1-19 GeV [l^. We use also following values 
of light quark masses: md(Ax) = 10 MeV, m„(Ax) = 7 MeV jT^ . 

Eq. (|l^ contains sizeable theoretical uncertainities due to non-perturbative strong in- 



teractions. However, as we show in the next section, for most parameter choices E^ alone 
gives the leading contribution to the neutron EDM. Therefore, one may hope that those un- 
certainities affect mainly the overall normalization of the neutron EDM. They do not affect 
significantly the possible cancellations between the phases (or in their coefficients), as long 
as such cancellations must occur predominantly inside the E^. For instance, if one uses rela- 
tivistic quark-parton model ||15|, |] instead of the naive chiral quark model, the contribution 
of quarks to the neutron EDM is weighted by |jl5|] their contributions Ag to the spin of pro- 
ton, measured in polarized lepton-nucleon scattering: En = rje^AuEd + A^E^ + AgEs). For 
the neutron, applying the isospin relations, we have (A„)„ = A^^, {Ad)n = A^, (As)„ = A^. 
With the assumption of degenerate masses and phases of the A parameters of the first two 
generations of squarks, Ed and Es have exactly the same structure and differ only by the 
quark mass multiplying the appropriate loop integrals: Ed/rrid = Es/mg (see eq. ( |B.30|) ). 
Therefore, accepting the values given in ref. [|l^: A„ = 0.82, A^ = —0.44, A^ = —0.11, one 
can estimate that, in the quark-parton model, the contribution from the s sea-quark, not 
present at all in the naive chiral model, is predicted to dominate. The leading contributions 
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to the neutron EDM in the considered models are (we assume rris ^ 150 MeV): 

4 
Chiral quark model : En ~ izVeEd ~ l-SSr^e-Ed 

Tfl 

Quark parton model : En ~ (A^ H -/S.s)rjeEd ~ — O.SSr^e-E'd (15) 

rrid 

and the u quark contribution is smaller. The coefficients of Ed for the two models differ in 
magnitude and even in sign but in both models the limits on the phases can be avoided if the 
Ed itself vanishes due to some internal cancellations. Of course, one should remember that 
the estimate (0) is subject to some uncertainity in A^ and a larger cancellation between d 
and s quark contributions, leading (at least in the low tan/5 region [|15|) to the dominance 
of the u quark contribution, is not excluded. Thus, there are uncertainities related to the 
choice of a particular quark model and uncertainities within the chosen model (like e.g. 
Imy4(i = \m.As assumption). Any theoretical calculation of the overall normalization of the 
neutron EDM should be considered as a qualitative one. The pattern of cancellations we 
discuss in the naive chiral quark model can be more generally trusted to the extend the 
smallness of the u quark contribution compared to the d and s quark contributions can be 
trusted as a model independent feature. At present the limits on the phases given by the 
electron EDM are more precise and better established. 

It was recently pointed out [l^ that 2-loop contributions to the neutron EDM may 
be numerically significant, especially for large tan/5 regime. Unlike most of the terms in 
eq. ([l^ , they depend mainly on the masses and mixing parameters of the third generation 
of squarks. Therefore, they are especially important in the case of the third generation of 
squarks significantly lighter than the first two generation, so that the 1-loop contributions 
are suppressed. We do not include such corrections in the present analysis. 

The formulae for the up- and down quark electric dipole moment are the following: 

E'd = £-2J:i:rncM{iKuc)L\V,^^y^'^*)(Cn{ml^,mlJ + \^^^^ 
°^ j=ik=i ^ "^ 



4: U 



—^\M,\ E lm{Zh',ZJ,'^)Cu{ml^, mi') (16) 

^^ fc=i 



K = -^EE^c,Im((V;^^)f^(V;^^)r)(Cn(m^^,m|j + iCi2(m|^,m^; 

g 4 6 . 

+ T^J-^ E E "^ivolm ( {V^uN)L\KuN)R^*)Ci2{ml^,m%o 
^^^ j=ik=i ^ ^ 
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+ ^l^3lEM41^^'K*)C^i2(<,|M3p) (17) 

^''^ k=l 

Chromoelectric dipole moments of the quarks are given by: 

-LO^ j = lk=l 
4 6 

+ T7h E E "^ivolm {{VaDN)L\KDN)R^*) Ci2(m^^, m%o) 
J-"^ j=lfc=l ■' 






^.' = T|r^EE"^c,Im((K^^)f^(K^^)^^^-^)Ci2K^,m^; 

-LO''' j = \k=\ 
4 6 

-LO''' j = lfc=l ^ 

"s I Jl ^' I V^ T,^ /' rylk rylk*\ o /^ /'IS/flS ,^„2 \ i 



Msl E Im(414'i?) (3Cn(|M3p,m|J + -C^^K^, |M3n) (19) 



i^' .=1 



Finally, the gluon chromoelectric dipole moment is given by: 

3a'^a^m^ ,« ,« m^ — mj ( mj mj rnl \ 

C, ^ ^|^MZ-Z»«)^^ffy|,.^.^ (20) 



where the definition of the 2- loop function H can be found in |T^ . 

The full form of all necessary fermion-sfermion-chargino/neutralino vertices can be found 
Appendix A| . In [Appendix B| we give also mass insertion expressions for the quark electric 



m 



and chromoelectric dipole moments. Although somewhat more complicated than in the case 
of leptons, they are very useful for qualitative understanding of the neutron EDM behaviour. 

3.2 Limits on phases 

The neutron EDM depends on more phases than the electron EDM. All electric and chromo- 
electric dipole moments depend on the common yU phase, but some of them are proportional 
to /i tan/? and others to ^ cot /?, hence the limit on /i phase does not scale simply like 1/ tan/5, 
as it was in the electron case. In addition, the quark moments depend on the phases of the 
two LR mixing parameters of the first generation of squarks. Ad and Au- The gluonic chro- 
moelectric dipole moment depends in principle on all A parameters and squark masses, but 
contributions from different squark generation are weighted by the respective fermion mass, 
so we take into account only the dominant stop contribution, dependent on the At parameter. 
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In practice, the analysis of the dependence of the neutron EDM on SUSY parameters 
appears less complicated than suggested by the above list, as some of the parameters have 
small numerical importance. As discussed below, the result is most sensitive to the squark 
masses of the first generation (left and right), gaugino masses, |/i|, A^, A^ and tan/3. 

The number of free parameters can be further reduced by assuming GUT unification with 
universal boundary conditions. Such a variant was thoroughly discussed in 0, §], so we do 
not repeat the full RGE analysis here, however its results can be qualitatively read also 
from the figures presented in this section. This can be done with the use of the following 
observations: 

i) As mentioned above, the neutron EDM is sensitive mostly to the masses of the first 
generation of squarks. Assuming universal sfermion masses at the GUT scale one 
can to a good approximation keep them degenerate also at Mz scale. The remnant 



m 



u 



of the GUT evolution is their relation to the gaugino masses: rriQ ^ mjj 

tUq + 6.5M^/2 ~ rriQ + 10|M2p, which leads to the relation rriQ ^ mjj k, rriD > 3M2. 

ii) The /i phase does not run. The /i itself runs weakly, it means that Im/i also runs 
weakly. It is a free parameter anyway. 

in) The imaginary parts of the first generation A parameters, ImA^ and ImA^, do not 
run, apart from the small corrections proportional to the Yukawa couplings of light 
fermions. Real parts of Au and Ad run approximately in the same way. Therefore 
universal boundary conditions at the GUT scale lead simply to 0a„ = <pAd at the Mz 
scale. 

iv) RGE running suppresses the At phase (present in the chromoelectric dipole moment 
of gluons Cg). Therefore, the low energy constraints are easy to satisfy even with large 
(f)At at the GUT scale. The limits on (f)At at the electroweak scale appear themeselves 
to be rather weak. 

v) With universal gaugino masses and phases, Mi/ai = M2/a2 = M^/a^, the common 
gaugino phase can be completely rotated away. 

Using i)-v) one can use our plots for the universal GUT case, just assuming common A 
phase, neglecting (pAt and looking at the part of plots for which rriQ > 3M2. Actually, in 
all Figures of this Section we keep degenerate squark mass parameters Mq = M£, = Mjj, so 
that the physical masses differ by D-terms only. We plot the results in terms of the physical 
mass of the Z)-squark rriD. 
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We consider first the limits on the /i phase, neglecting the possibility oi fi—A cancellations. 
In Fig. |16| (analogous to Fig. ^) we show where the generic limits for the fi phase given by 
the neutron EDM are strong. We plot there the area where the limit on | sin0^| tan/? given 
separately by each of the contributions present in eq. QT^) is stronger then 0.2 or 0.05 (in 
addition we consider separately the chargino, neutralino and gluino contributions to Ed). 
For small |/i|, IM2I, squark masses mq ~ mo ~ ttlu > 1600(750) GeV are required to avoid 
the assumed limits. 
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Figure 16: Regions for which generic limits on |sin0^|tan/? given by neutron EDM are 
stronger then, respecively, 0.2 (left plot) and 0.05 (right plot). Degenerate squark masses 
and GUT related Mi, M2, M3 assumed. 



The dominant contributions to the coefficient multiplying sin 0^ come from the first term 
of eq. (|14D) i-6. from the d-quark electric dipole moment, as illustrated in Fig. |l^. The only 
exception is large |/i| and light gauginos case, where also Cd becomes comparable to the 
other term. In principle, relative importance of various contributions changes with tan/?, 
as Ed, Cd are proportional to /itan/3 and E^, C^, Cg to ji cot (3. However, Ed and eventually 
Cd always dominate and the /x phase coefficient scales again, like in the electron case, as 
tan/3 (Fig. |1^ has been done for tan/? = 2). The largest contributions to ^ phase coefficient 
are given by the chargino and gluino (for small and large |/i|, respectively) diagrams. They 
have the same sign, so again, like in the electron case, the total /i phase coefficient may 
disappear only if one allows the non-universal gaugino phases. The limits on | sin0^| tan/5 



on niD — 1^2 1 plane are plotted in Fig. 18 



Some differences with the electron case may be observed in the structure of possible 
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Figure 17: Contributions to sin0^tan/9 coefficient in the neutron EDM, normalized to (di- 
vided by) the experimental limit. Solid, dashed, dotted line: chargino, neutralino and gluino 
contributions to Ed] dotted-dashed line: E^ + Cd + Cu + Cg summed up. IM2I = 100 GeV, 
Mi,M2,M3 GUT-related, 0^. = and degenerate squark masses assumed. 
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and degenerate squark masses assumed. 
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H — A cancellations. For E^ the term proportional to A^ originates from the neutralino 
exchange diagram. For the neutron, additional contributions proportional to A^, A^ and 
At are given by the diagrams with gluino exchange and they have larger magnitude than 
those induced by neutralino loops, as illustrated in Fig. |19| (this efect is particularly strong 
for large |/i| and light gauginos). This means that, on the one hand constraints on the Aj 
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Figure 19: Coefficients of sin 0^ tan /5, \Ad\/m]j sin (j)Aa, |A„|/m(/sin0^^ and \ At\ /rriT sin (p At 
terms in the neutron EDM (solid, dashed, dotted and dashed-dotted lines respectively). 
Ml, M2, M3 GUT-related and degenerate squark masses assumed. 



phases are somewhat stronger than in the electron case but, on the other hand, smaller Aj 
values are necessary for cancellations. For small |/i| ~ 100 GeV one needs Ae/rriE > 14 but 
only Ad/rriD ~ Au/mu > 3. Furthermore, in unconstrained MSSM we have bigger freedom 
because of several different Ai parameters present in the formulae for En- Contributions 
proportional to A^jAt, coming from E^, Cu and Cg, are more important comparing to those 
given by the E^, Cd as they are not suppressed by the relative factor cot^ /5 (like it is for 
the fi phase). Therefore, one has to take into account all A^ phases. In Fig. ^ (compare 
Fig. ^) we plot the regions of mD — \Mi\ plane allowed by the neutron EDM measurement 
assuming 0^ = 0^^ = 0a„ = 7r/2 and various values of Au/mu = Ad/mu. 

The overall conlusion is that eventual cancellations in neutron EDM are more likely than 
in the electron case. They require somewhat smaller values of A parameters when one 
considers fi — A phases cancellations. Furthermore, assuming non-universal A^ parameters 
it is possible to suppress simultaneously both Ee and En values below the experimental 
constraints, at the cost of rather strong fine-tuning if the SUSY mass parameters are light. 
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4 /J, phase dependence of ex , Am^ and b 
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Analysing the dependence of K^K'^ and BB mixing on the SUSY phases, we assume that 
there is no flavour violation in the squark mass matrices, so that only chargino and charged 
Higgs contributions to the matrix element do not vanish (gluino and neutralino contributions 
are always proportional to the flavour off-diagonal entries in the squark mass matrices). 
Furthermore, only chargino exchange contribution depends on the yU, M2 and A phases and 
is interesting for our analysis. The leading chargino contribution is proportional to the 
d^^^^dj^ ■ dif'-ffj,dj^ matrix element and has the form: 



1 



{Mc)llll = - E iV,ocynV,ocyL\V,oc)f'iV,ocYL'''D2iml^, m\,m\,ml) (21) 

where one should put / = 2, J = 1 for K^K^ mixing, / = 3, J = 1 for B^Ba mixing and 
/ = 3, J = 2 for BgBs mixing (see [Appendix A] for the expression for loop function D2). 



In order to analyze the dependence of the matrix element ( pT]) on the phases, we consider 
the simplest case of flavour-diagonal and degenerate up-squark mass and L-R mixing matrices 
and |/i|, IM2I > 2Mz. In this case we can expand the matrix element in the mass insertion 
approximation, both in the sfermion and chargino sectors, as described in [Appendix B| . The 
eq. (PH) gives in such approximation: 



(Mc) 



LLLL 



8 



-{K^Y^K)^^[D2{\^\\\fi 



mjj,mjj} 
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+ 8M^Re [(/i* cos 13 + M2 sin {3) {fx cos p + Afj sin /?)] 



X 



5 ^2(|/i|M/i| 



,mu,mu, 



-D,{\fx\MM,\^ml,ml) 



drnfj 



|M,I2 



(22) 

e/^ and Atub are proportional, respectively, to the imaginary and real part of the matrix 
element. One can see immediately from the equation above that in the leading order it is 
sensitive only to \fi\ and to the real parts of the M2/i, Aufi and M2A\j products, i.e. to 
cosins of the appropriate phase combinations, not sins of them like it was in the EDM case. 
Eventual effects of the phases can be thus visible only for large phase values. Even then, 
they are suppressed by M^/mf^ ratio and small numerical coefficient mutiplying them. An 



example of the ex dependence on the /x and Au phases (based on exact expression (^) and 
assuming M2 to be real) is presented in Fig. ^ As can be seen from the Figure, even for 
light SUSY particle masses the change of the ex value with variation of /i and A phases is 
smaller than 5%. 
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Recently the potential importance of non- leading chargino contributions to ex and BB mi- 
xing (proportional to LLRR and RRRR matrix elements) was discussed in the literature |T^ . 
Such contributions are suppressed by the factor m^(mf) tan^ (3/My^ for K^K^ {BB ) mixing. 
Therefore they can be significant only for large tan j3 values. Demir et al. [|l^ estimate that 
ex ~ 3 ■ 10^^ can be obtained solely due to the supersymmetric phases of fi and At, as- 
suming vanishing Kobayashi-Maskawa phase 6xm = 0. It requires however large tan/3 = 60, 
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6^j ~ —7t/2 and hght stop and chargino: in the example they give 



l/^l 



I Mo 



125 GeV, m 



tL 



mx 



150 GeV, |Ai| = 250 GeV, what gives physical masses 



mx 



83 GeV, r/ici ~ 80 GeV. As follows from our discussion in sections ^ and |^, it is rather 
unhkely, although not completely impossible, to avoid limits on /i phase for light SUSY spec- 
trum and simultaneously large tan/3 (limits on /i phase are inversely proportional to tan/9). 
The generic limits are in such a case very tight, so in order to avoid them one needs strong 
fine tuning and large cancellations between various contributions. For example, assuming the 
above set of parameters, universal gaugino masses and phases, uie = mu = rrijj = 200 GeV, 
Au = Ad, 4>A, = 7r/2 one needs very large \Ae\/mE ~ 750 and \Au\/mu = \Ad\/m£, ~ 170 in 
order to avoid the limits from the electron and neutron EDM measurements. Alternatively, 
one can keep smaller Ai but give up the assumption of gaugino mass universality and adjust 
precisely their masses and relative phases. In each case, a very peculiar choice of parame- 
ters is required to fulfill simultaneously experimental measurements of e^- and EDM's if one 
choose to keep 5km = 00. 

Because of the weak dependence of ex and Atub on flavour conserving supersymmetric 
CP phases (excluding possibly very large tan j3 case) , the determination of the KM matrix 



^Bounds on the /i phase cannot be avoided even by assuming heavy first two generations of sfermions. 
In this case, for large tan/3 considered here, 2-loop contributions from the third generation discussed in 
ref. jl^ are themselves sufficient to exclude the large /i phase necessary to reproduce the result for €k - see 
the discussion in lEOl. 
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phase 5km (see e.g. [^) is basically unaffected by their eventual presence. However, one 
should remeber that the KM phase detemination in the MSSM depends on the chargino and 
charged Higgs masses and mixing angles, even if not on their phases. 

In contrast to e^ and Aitlb, b —>■ sj decay appears to depend strongly on the fi and At 
phases. This is illustrated in Fig. ^ Exact formulae for the matrix element for b ^ s'y decay 
can be found e.g. in ||21|. Expansion in the sfermion and chargino/neutralino mass insertions 
(assuming degenerate stop masses) gives the following approximate result for the so-called 
Cj coefficient, which primarily defines the size of Br(B —^ Xsj): 

rio)(M ^ ^ 1 f ("^h ^^ f'( ""^ ^ M'' cot 13 + At^x) m} 

12|/ipsm p |/i|^ o|iW2p 1^2! o|/ipsm/p |/ip 

M^Re[(/icot/5 + A^)(/i^cot/? + M2)] [ /((p#) f'Aw 



3(|M2|2-|/i 



M'w 



3(|M2|2-|/i|2) 



fl( "t \ f'(—L\ \ / f'f i \ f'( 

/2I1M2PJ /2l|^|2j I ig^-^p^^M^ ' ^^l^2P^ •^2l|^| 



IM2I2 |/i|2 •- ' ^ ^\ |j^^^|4 |^|4 




(23) 
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where 



, , , x{?,x - 2) , 8x2 + 5x - 7 

(1 — xp (1 — xY 

—x{?>x — 2) , x(5x — 3) 

; -;— log X -, 

'l-xY ^ 2 1-x 



/2(^) = ;; .3 l^g^-^7?— ^ (24) 



Expression (^) contains terms proportional to both real and imaginary parts of the /i and 
At parameters. However, the branching ratio Br(B -^ Xs'j) depends on \C^ (Mw)]"^ which, 
like in the ex case, depends mainly on the real parts of the /i and At parameters, i.e. on 
cosins of the phases. This is clearly visible in Fig. |2^. However, contrary to the ex case, this 
dependence is quite strong and growing with increase of tan/3 and of the stop LR- mixing 
At parameter. Also, as follows from the discussion in the previous Section, the limits on At 
phase are rather weak, independently on tan (3, so one can expect large effects of this phase 
in Br{B — > ^^7) decay. 

5 Conclusions 

In this paper we have reanalyzed the constraints on the phases of fiavour conserving super- 
symmetric couplings that follow from the electron and neutron EDM measurements. Also, 
we have discussed the dependence on those phases of ex , Ams and the branching ratio for 
b ^ sj . We find that the constraints on the phases (particularly on the phase of fi and of 
the gaugino masses) are generically strong < 10^^ if all relevant supersymmetric masses 
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are light, say < (9(500 GeV). However, we also find that the constraints disappear or are 
substantially relaxed if just one of those masses, e.g. slepton mass, is large, tjie > 0{1 TeV). 
Thus, the phases can be large even if some masses, e.g. the chargino masses, are small. 

In the parameter range where the constraints are generically strong, there exist fine-tuned 
regions where cancellations between different contributions to the EDM can occur even for 
large phases. However, such cancellations have no obvious underlying symmetry principle. 
From the low energy point of view they look purely accidental and require not only fi — A, 
fj, — Mgaugino OT Mi — M2 phasc adjustmcut but also strongly correlated with the phases 
and among themselves values of soft mass parameters. Therefore, with all soft masses, say, 
< (9(1 TeV) models with small phases look like the easiest solution to the experimental 
EDM constraints. This conclusion becomes stronger the higher is the value of tan (3, as the 
constraints on fi phase scale as l/tan/3, and will be substantially stronger also for low tan/3 
after order of magnitude improvements in the experimental limits on EDM's. Nevertheless, 
since the notion of fine tuning is not precise, particularly from the point of view of GUT 
models, it is not totally inconceivable that the rationale for large cancellations exists in the 
large energy scale physics (in a very recent paper [^ it is pointed out that non-universal 
gaugino phases necessary (but not sufficient) for large cancellations can be obtained in some 
String I type models). Therefore all experimental bounds on the super symmetric parameters. 



and particularly on the Higgs boson masses |23], should include the possibility of large phases 



even if with large cancellations, to claim full model independence. 

The dependence of ex and AniB on the supersymmetric phases is weak and gives no 
clue about their values. Hence, the 6km determination remains essentially unaffected. 
Large effects may be observed in b —>■ s'j decay, but, apart from the 0^ and (pAt phases, 
b ^ s'y amplitude depends on many free mass parameters, including the Higgs mass and the 
masses of squarks of the third generation. 

Appendix A Conventions and Feynman rules 

For easy comparison with other references we spell out our conventions. They are similar to 



the ones used in ref. |23- The MSSM matter fields form chiral left-handed superfields in the 
following representations of the SU{3) x SU{2) x f/(l) gauge group (the generation index is 
suppressed) : 
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Scalar field Weyl Fermion field SU{3) x SU{2) x f/(l) representation 
L^(%) l^{:) (0.2,-1, 

E^ e" (0,0,2) 

D^ d" (3,0, 1) 

[/'= «= (3,0,-|) 

Two S'f/(2)-doublets can be contracted into an S'f/(2)-singlet, e.g. H^H^ = eijHjH'^ = 
—H\H2 + H2HI (we choose ei2 = —1; lower indices (when present) will label components 
of S'?7(2)-doublets). The superpotential and the soft terms are defined as: 

W = Wq + Wcp (A.l) 

Wo = Y^H^LE + YdH^QD + YuH'^QU (A.2) 

Wcp = fiH^H^ (A.3) 

^soft = ^soft~0 + ^soft-CP (A. 4) 

-Q^MlQ - D^^MlD^ - W^MIjW (A.5) 

+YeA,H'LE' + r.Ad/jiQ/}^ + Y^A^H^QU' + H.c. (A.6) 

where we divided all terms into two subclasses, collecting in Wcp and Cgoft-cp those of 
them which may contain flavour-diagonal CP breaking phases. We also extracted Yukawa 
coupling matrices from the definition of the Aj coefficients. 

In general, the Yukawa couplings and the masses are matrices in the flavour space. Ro- 
tating the fermion fields one can diagonalize the Yukawa couplings (and simultaneously 
fermion mass matrices). This procedure is well known from the Standard Model: it leads 
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to the appearance of the Kobayashi-Maskawa matrix K in the charged current vertices. In 
the MSSM, simultaneous parallel rotations of the fermion and sfermion fields from the same 
supermultiplets lead to so-called "super-KM" basis, with flavour diagonal Yukawa couplings 
and neutral current fermion and sfermion vertices. As we neglect flavour violation effects in 



this paper, we give all the expressions already in the super-KM basis (see e.g. |21] for a more 
detailed discussion). 

The diagonal fermion mass matrices are connected with Yukawa matrices by the formulae: 



mi 






m,. 






rrid 
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The physical Dirac chargino and Majorana neutralino eigenstates are linear combinations 
of left-handed Winos, Binos and Higgsinos 



a 



iZ!}W- + Z'Jhf 



(A.7) 



where W^ = {W^ T W^)/V2. 



N' 



-iZj^*B - iZ'^*W^ + Z'^*h\ + Z'^*hl 



(A.8) 



V iZ'^B + iZ'^W^ + Z'^h\ + Z'^hl J 
The unitary transformations Z^, Z^ and Zn diagonalize the mass matrices of these fields 
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(A.9) 



and 
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4-component gluino field g is defined as 





M2 
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2 



2 
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^ -^^ 



g'v2 \ 
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-/i 




^N 



(A.IO) 
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(A.ll) 



The sfermion mass matrices in the super-KM basis have the following form: 



Ml 



cos 2/3 , 

mu(cot/?/i*l + Au) 



KMlK^ + ml-^^{M_ 



4M2,)1 



-m. 



„(cot/3/il + A^) 



Ml + ml + ^^{Ml - M^)l 



Ml 



-md(tan/3/i*l + Ad) 



— m(i(tan/5/il + A*^ 
Ml + mj-^{M] 



a 



M^)l 
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^2 _ I M| + m2 + ^(M|-2M2,)l -m,(tan/?/il + A^ 



'L 



-m, 



,{tanpi^n + Ae) Ml + ml- cos2/5(M| - M^,)! 



Ml = Ml + ^M|i (A.12) 

where 6w is the Weinberg angle and 1 stands for the 3x3 unit matrix. 

Throughout this paper we assume that there is no flavour and CP violation due to the 
flavour mixing in the sfermion mass matrices, i.e. matrices Mf, Aj are diagonal in the super- 
KM basis. However, one should remember that in this basis mass matrices of the left up and 
down squarks are connected due to gauge invariance [^. This means that it is impossible 
to set all the {M'Jj)ll to zero simultaneously, unless Mq ~ 1. 

The matrices M.1, -^i) -^^ ^ind M.'^j^ can be diagonalized by additional unitary matrices 
Zy (3 X 3) and Z^, Zjj^ Zd (6 x 6), respectively 



Mi) =ZlMiZL (A<y =Z^j,M'f,Zj, (A.13) 

The physical (mass eigenstates) sfermions are then defined in terms of super-KM basis 
fields (^) as: 

v = Z% L = Zl(^A U = z\j(llA D = Z^^(I^A (A.14) 

In order to compactify notation, it is convenient to split 6x6 matrices Z^, Zjj, Z^ into 
3x6 sub-blocks: 

Z^^Z^^ 4^''^ = ^XR /=1..3,^ = 1..6 (A.15) 

(the index i numbers the physical states). Formally Zxl and Zxr are projecting respectively 
left and right sfermion fields in the super-KM basis into mass eigenstate fields. 

Using the notation of this Appendix, one gets for the Feynman rules in the mass eigenstate 
basis: 

1) Charged current vertices: 
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V' 



d' 



W 



w" 



D: 



{^C+Y -i[g2Z+PL + Yi^Z^*PR 



7lJ* 



{ctr ' 



-g,Z(jlZt,+Y;^Zi^*^ZtAPL 






K 



ji 



-CI 



- {g^ZJ^lZt, + YjZJ^^^Z^;) Pl + Y,'Zj^lZ+;Pn] K^^ 



2) Neutral current vertices: 



L- 



u 



ut- 



■iVO 



-iV° 



'^^LL [9iZn + 92Zn) + ^i Zj^^Zj^j Pl 

+ [-grV^ZtRZ]^"^ + Y.^Z^lZT) Pr 
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d' 



D- 



{^,Z^d1 (f ^^^' - 92Z'^) + YlZ^^^Z"^) P, 
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?1V^ T'/j* 7^3 



^DR^N + ^d ^DL^N ) ^R 



3) Neutral current vertices with strong coupling constant: 
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U- 



-tgsV2T:^ {zIJIPl - ZIj%Pr 
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, /OT'1 ( yli* P yli* p 
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Finally, we give here explicit formulae for the loop functions. Three point functions 
Cii,Ci2 are defined as: 



p,g-*0 



d^k 




k^" 




(27r)4 [A;2 - 


- mf][(p + ky - 


- mf\[{k + p + qY - 


- 7712] 
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{Any 



(p^CuimlmD+q'^Cuimlml)) (A.16) 



A[x — y)"^ 2[x — y)'^ x 
n f \ ^^y xy y 

The four point loop function D2 has the form 

n f \ ^^ 1 "^ y'^ ■[ y 

{x — v){x — y){x — z) V {y — v){y — x){y — z) v 



^2 



log- (A. 19) 



{z — v){z — x) {z — y) V 

Appendix B Feynman rules for mass insertion calcu- 
lations 

We list below the Feynman rules necessary to calculate contributions to lepton EDM in the 
mass insertion approximation. We treat now off-diagonal terms in chargino and neutralino 
mass matrices (proportional to ^1,^2) as the mass insertions. In order to reduce remaining 
charged and neutral SUSY fermion mass terms to their canonical forms: 

- \fi\h^h- - h^^\Wh°^ - \M2\W^W~ - ^\M2\W^W^ - huilB^B^ (B.20) 

we define the 4-component spinor fields in terms of the initial 2-component spinors as: 

^' = ^[ \-J.-^A 5° = M -J.-^\ (B.22) 



-e"*' 



</>2/2p^3 / I _g-*0i/2^ 






Then the necessary Feynman rules are: 

, E{E'') 
p i ±ie{p + ky 

I 

7 k E{E^) 
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h^ 
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/lO 1^0,^0 



i^i _!,' \{v,-v,){e-^('^^+'^^^-y'PL 
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_gi{0M+0i,2)/2p^'j 
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-^f- 
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^ (t^i^L + t;2e^(<^''+'^^)Pi? 



X 



-^{ViAe + V2fi*Ye] 



E E^ 

As a cross check of a correctness of the calculations and a comparison of exact result with 
direct calculation of diagrams with mass insertions one can also expand the exact results in 
the regime of almost degenerate sfermion masses and large fi and Mi^2- 

For sfermion masses we assume that they are almost degenerate and expand them around 
the central value (X = u, L, D, U): 



"^i 



m^ + 5m j^ 



We expand also the loop integrals depending on the sfermion masses: 



f{m\_,m\,...) ^ f{m\,m\,...) + {m\,-m\) 



dm% 
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X, 
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+ ... 



[B.2A) 



df{m\,m\,...) 



m - =m ■ 



(B.25) 



and use the definitions of the mixing matrices Zx and their unitarity to simplify expression 
containing combinations of the sfermion mixing angles: 



Z'^Zj^* = 6'^ 



Z^Z^^m^^ 



Ml 



;B.26) 



In order to see the direct (although approximate) dependence of the matrix elements with 
chargino exchanges involved on the input parameters such as /z and M2, one should expand 
also chargino masses and mixing matrices. Such an expansion can be done assuming that 
M2, /i and also their difference are of the order of some scale |/i|, IM2I, |/U — M2I ~ A ^ Mz- 
Actually corrections to the most physical quantities start from the order (!?(M|/A^), so 
already A ~ 2Mz leads to reasonably good approximation. 
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Matrices Z and Z"*", diagonalizing the chargino mass matrix (eq. 



may be chosen 



as 



1 -a* 
a 1 



-l<p2 



e 



-(3*6 



-icpfj 



-i4>ii 



where 



a 



e V2M2 + f i/i* 



/3 



e V1M2 + V2^* 



(B.27) 



(B.28) 



V2sw IM2P - l/iP ^ y2svv^ IM2I2 - |/i|2 

(in the above expressions we do not assume M2 to be real). Phases in Z^, Z^ are chosen to 
keep physical chargino masses nic^, mcj in eq. (|A.7| ) real and positive: physical results do 
not depend on this particular choice. 

Eqs. ( |A.7| , p.27|) give IM2I and |/i| as the approximate masses for charginos: 
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;b.29) 



Expansion of the neutralino mass matrix is more complicated and tricky as two Higgsinos 
are in the lowest order degenerate in mass. The appropriate expressions can be found in |2^ . 
In this case direct mass insertion calculation using the Feynman rules given in this Appendix 
is usually easier. 

Below we list mass insertion approximation expressions for the electric and chromoelectric 
dipole moments of d and u quarks. We assume equal masses of the left and right sfermion of 
the first generation mu ^ muc ^ m^i ^ tti/jc = mq and neglect small terms proportional to 
higher powers of the Yukawa couplings of light quarks. In such a case, corresponding dipole 
moments can be approximately written down as: 
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+ 3(4;0^I-(^^^)^°^^ W-\M,\- 

^(i;^I"^[^3(^cot/5 + AJ] ^-2 (B.31) 

Analogously, chromoelectric dipole moments of quarks read approximately as: 



^'^ " ^(4^Im(M2/.)tan^ H^^P^P 

+ 



2(47r)2 ^ ^^' ' |^|2-|Mi|2 
^^^-^lm[M,(Man^ + A,)] ^-^ 

^^^Im[M3(Man^ + A,)] (^3 ^ + - ^^ j (B.32) 



3^2 „ Ci2(m|,|/in-Ci2(m^,|M2p) 

^" " ^(i^I-(^^/^)-t/^ |^p_|M2P 

^ /^^m^, .,. . ,^ Ci2(m^,|/ip)-Ci2(m^,|Mip) 

+ ^(i;O^I-(^^^)^°^^ |;.p-|M,P 

2^'^^-^-Im[M,(,cot^ + A:)]^^-«'l^^l^) 



9(47r)2 ^ '^^ /- ' «/j ^9 



|2^ 



+ ^^I-[^3(/.cot/3 + AJ](^3 ^^^^j + - ^^ J(B.33) 
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